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Proposals for Bell inequality tests on systems restricted by the particle number superselection 
rule often require operations that are difficult to implement in practice. In this paper, we derive 
a new Bell inequality, where measurements on pairs of states are used as a method to by-pass 
this superselection rule. In particular, we focus on mode entanglement of an arbitrary number of 
massive particles and show that our Bell inequality detects the entanglement in an identical pair of 
states when other inequalities fail. However, as the number of particles in the system increases, the 
violation of our Bell inequality decreases due to the restriction in the measurement space caused by 
the superselection rule. This Bell test can be implemented using techniques that are routinely used 
in current experiments. 

PACS numbers: 



I. INTRODUCTION 

The rise of quantum information science over the last 
decade has made it popular to seek and understand quan- 
tum many-body systems that contain entanglement [T]. 
One group of these many-body systems are the ultra- 
cold atomic gases, such as Bose Einstein condensates 
However, unlike spin lattice systems where the par- 
ticles are distinguishable through position, particles in 
ultra-cold gases are typically indistinguishable from one 
another. Indistinguishability means that the first quan- 
tised many-body state of the particles should be (anti)- 
symmetrized, but then the corresponding state space no 
longer has the tensor product structure required to define 
entanglement [3]. 

One can, however, recover a tensor product state space 
by viewing the ultra-cold gases in terms of second quan- 
tised modes [1]. It has been shown that entanglement 
naturally exists between spatial modes in non-interacting 
Bose Einstein condensates [SHZj and in other ultra-cold 
atomic gases provided the coherence length of the par- 
ticles extends over the selected modes f5]. In order to 
perform a Bell test on such systems, the spatial modes, 
which behave in some sense like a pair of qudits, must 
be rotated away from the particle number basis. How- 
ever, since systems of massive particles are restricted by 
a superselection rule [I7HI9] that forbids rotations away 
from the subspace of fixed particle number, a Bell-like 
test of the mode entanglement of massive particles is not 
straightforward. On the other hand, spatial mode entan- 
glement (and non-locality) of a single photon has been 
extensively studied and the experimental verifica- 

tion of single photon entanglement has been obtained [T5] 
via the CHSH Bell test 

Despite the superselection rule, a few schemes to test 
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a Bell inequality with a single massive particle have re- 
cently been put forward In [H] it was suggested 
that the spatial modes could be rotated away from the 
particle number basis by coupling to a coherent particle 
reservoir. However, to reliably confirm spatial mode en- 
tanglement of a single massive particle in an experiment, 
one would have to ensure that no additional entangle- 
ment entered the system via the particle reservoir (see 
P51 [23] for further discussion of this point) and even if 
this were guaranteed, such a coupling to the reservoir 
is likely to be difficult to implement in realistic condi- 
tions. Moreover, in both the photon experiment |15j and 
in theoretical proposals [HI HI] , a local post-selection of 
the measurement outcomes was required to see a vio- 
lation of the CHSH inequality, which meant that these 
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FIG. 1: Basic set up for the Bell test of mode entanglement 
using two copies, p^'^' ® of a bi-mode state with A*' 

and A^' massive particles respectively. Modes a and A are 
given to Alice and 6 and B to Bob. Each party makes a 
general measurement in the subspace of fixed particle number 
by sending their modes through a beamsplitter parameterised 
by the local angles (pA for Alice and tps for Bob. 
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tests only strictly probed for entanglement (i.e. the Bell 
inequality was used as an entanglement witness [25j). and 
not the more stringent property of non-locality. 

In this paper, we derive a general Bell inequality to 
detect the mode entanglement of massive bosons, which 
overcomes the possible pitfalls pointed out above. In par- 
ticular, we use two copies of the system to make rotated 
measurements despite the superselection rule. Using two 
copies not only eliminates the need to have careful cou- 
plings to a particle reservoir in order to overcome the 
superselection rule, but also ensures that no auxiliary 
entanglement from this reservoir will be responsible for 
the violation. Note that Wiseman and Vacarro consid- 
ered the entropic properties of pairs of superselection rule 
restricted states in [T^ . Unlike previous tests of spatial 
mode entanglement that focus solely on single particles, 
our Bell inequality can be applied to systems contain- 
ing any number of bosons, which can be both massless 
or massive. At present there has been no experimental 
test for the existence of mode entanglement of massive 
particles - even though this type of entanglement is pre- 
dicted to be ubiquitous in ultra-cold atomic gases and 
that it has recently been shown to be useful for quantum 
communication 26J . Because our test is relatively simple 
compared to previous proposals to detect mode entangle- 
ment of massive particles, its implementation could allow 
for the first direct confirmation of mode entanglement of 
massive particles. 

We begin in the next section by explaining in more 
detail why two copies of a bi-mode quantum state allows 
one to make the measurements required to implement our 
Bell inequality despite the superselection rule. Building 
on this, we derive in section (III) a general Bell inequality 



for two copies of a bi-mode state each with N particles. 
In section (IV) we test our Bell inequality with some 



examples of mode entangled states and in section (|V|) we 
discuss these results. Finally, in section (VI) we suggest 



how one could implement the test in practice, referring 
to current experimental techniques that could be put to 



II. THE CHSH BELL INEQUALITY AND THE 
PARTICLE NUMBER SUPERSELECTION RULE 

Bell inequalities allow to test the remarkable ability for 
entangled states to violate local realism. Local realistic 
theories impose constraints on the correlations between 
measurement outcomes on two separated systems. For 
systems occupying the state space, H = ® C^, the 
most commonly tested Bell inequality [27] is the CHSH 
inequality [IS]. The CHSH inequality can be expressed 
in terms of joint expectation values of observables for two 
parties, A and B, as 

BCHSH = l(ilSi) + {A1B2) + {A2B1) - {A2B2)\ < 2, 

where {AiBj) = ti[AiBjPAB\ and Ai and Bj each have 
two outcomes for i, j = 1, 2. When this inequality is vio- 



lated there is entanglement between the two subsystems. 

In systems of massive particles described by quantum 
mechanics, particle number is a superselected quantity 
so that the total particle number operator, TV, com- 
mutes with all other physical observables. The system 
density operator, p, therefore cannot contain any off- 
diagonal terms that connect states of different particle 
number. The corresponding Hilbert space, "H, is decom- 
posed asH = (B'^^q'Hn, where Hn the subspace con- 
taining N particles. For systems of two spatially sep- 
arated modes, A and B, the state space Hn has the 

substructure Hn = ®,T=o('^n ® f^N-n)^ where Hn^^'' 
denotes the Hilbert space of mode A{B) with n parti- 
cles. 

Now consider a single copy of a two mode system con- 
taining one massive particle. A general state in the 
subspace. Hi, spanned by |01) and |10) can be writ- 
ten as pAB = p|V'"^)(V'"^Ub + ~ p)\ip~){i^'\AB, where 
\^^)ab = ^(|01)±|10)) and |01) = \0)a<E)\1)b denotes 
no particles in mode A and one particle in mode B. The 
parameter, p, determines the entanglement of pabj only 
when p = i is the state separable, with p = or 1 result- 
ing in a maximally entangled state. Because of the super- 
selection rule, the only valid local measurement is particle 
number and, if we were to test Eq. ([T]), every joint cor- 
relation function, {AiBj)p^^, would be (AiBj) p^^^ — — 1 
Since (AiBj) is independent of p, measurements 
in the particle number basis cannot distinguish separa- 
ble states from entangled ones, and such measurements 
never allow for violation of the Bell inequality, even with 
a mode entangled state. 

Conversely, we take now two copies of the above state, 
PxYi shared between two parties, a = Pab® Pcd^ where 
modes A and C are given to the first party and modes 
B and D to a second party. This composite state will 
allow to 'see' the entanglement of pxY despite the su- 
perselection rule. We consider now a toy example to 
highlight this point: Each observer will make an (incom- 
plete) local measurement on their two modes in the two- 
dimensional subspace of fixed particle number spanned 
by the normalised basis vectors, )}, where 

\ip+) = a|10) +/3e^'^|01) and |^-) = /3p>|10) - a|01). If 
the first party measures A{ip) ~ \ip-\-){ip + | — \ip—){ip — \ 
and the second party B{e) = \6+){d + | - \e-){e - |, 
the joint expectation values of Eq. ([T]) with the state, cr, 
have the form, (i(^)^(6l))^ = 8(p- \Y (3^ cos{lp - 9) . 
It is clear that the degree of correlation depends on the 
entanglement of the individual states, pxY, since the pa- 
rameter, p, is present. For instance, when the states 
are separable, i.e. p — the correlation function is 
{A{ip)I3{e))fj = 0, yet when p = and a = /3 = 

the correlation function is maximal, {A{(f)B{9))a = 
\ cos((/5 — 9). Moreover, the local parameters, and 6, 
can be altered by each party respectively to change be- 
tween different measurement settings (1 and 2) required 
for a Bell test. We will now expand this basic example to 
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formulate a general Bell inequality for pairs of N particle 
states. 



III. BELL INEQUALITY FOR PAIRS OF 
PARTICLE NUMBER SUPERSELECTION RULE 
RESTRICTED STATES 

Consider two systems split into two spatially non- 
overlapping modes, the first system has N massive parti- 
cles and the second N' massive particles. We will denote 
the total state of the system, with N + N' particles, by 
^(N+N') _ (g)p(^\ Party A (or Ahce) has access to 
modes, a and A, and party B (or Bob) to modes, b and 
B, see Fig ^. 

In our Bell test, Alice will make a joint measurement on 
her two modes in a subspace, Hm = ©*LoC^n«''^i-„): 
of fixed particle number, M, spanned by the basis states 

{|M, 0)aA, \M - 1, l)aA, • ■ • , |0, M)aA}, although WC do 

not know a priori what this number, Af (0 < M < 
A'' -t- N') will be. The operator for this arbitrary high- 
dimensional measurement basis is, for Alice, 

N+N' 

M^'a)^ X! e{nc,mc)\nc,mc){nc,mc\cc- (2) 

nc+mc—O 

Likewise, Bob will measure his two modes in the basis, 
Biip'i^), and will receive {N + N') — M particles. Here, 
e{nc,mc) are weighting coefficients. The local param- 
eters, ip^^ and (^^^ where i, j — 1, 2, denote the two 
measurement settings that Alice and Bob will use in the 
Bell test. Measuring each mode directly would allow for 
only local particle number measurements, but we would 
like to perform a general measurement within Hm- To 
do this each party passes their two spatial modes through 
a beamsplitter defined, for Alice, by the transformations 

c^aa + e"*^-^ fSA C ^ /3a - e-''^^ a A, (3) 

where c and C are annihilation operators for the two 
output ports. The operators a and A are annihilation 
operators for the two input modes of party A. There 
are similar beamsplitter transformations for Bob, where 
we denote the output modes as d and D. Each party 
measures the output modes in the particle number basis, 
the outcomes of which depend on the local parameter 
^PA(B)■ Hence, the number of particles, ric and mc, in 
the two output modes, c and (7, appear in the observable 
([2]) for Alice. A measurement of |nc,m(7)cC, corresponds 
to an effective measurement on the input modes, a and 
A, of 



(aat + (3e-"^^A^)""' 
\nc,mc)cc = ^=f X 



(4) 



(/3a 



a. 
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where |0,0)a a is the vacuum of modes a and A and we 
have used We will discuss in more detail a physical 
implementation of our test in section (VI). 

Since there are a total oi N + N' particles in the com- 
posite system, there are {1/2{N + N') + 1){N + N' + 1) 
number of different measurement outcomes. For exam- 
ple, when N = N' = 1 and for balanced beamsplitters, 
a = /? = ^ , we have a total of 6 outcomes: 

{nrnjcc measurement on the modes a and A e{nc,mc) 



1 00) 
1 10) 
101) 
1 20) 
111) 
|02) 



100) 



HOI)) 
HOI)) 



1(|20) -He-*^^A/2|11) -he-*2'^-*|02) 
^(|20)-e-«^-|02)) 



-1 



1(|20) - e'''fi^V2\ll) + e-'^'^^\02) 
The weighting function, e(nc, rnc), is explicitly chosen as 

(mc + ne)(,Tic+ne + l) 



e("c,"^c) = (-1) 



mc + - 



(5) 



which gives a sharp binning of results, as in |29| . There, 
sharp binning was shown not only to be optimal, but also 
to result in a tight Bell inequality. 

We construct a Bell-type inequality from the local ob- 
servables, A{ip^^^) and B{ip^^^). Like the CHSH type com- 
bination of Eq. (|l|, we formulate B — E{(p']l\ (p^j^'') + 
E{f^}\ V^B^ ) ^-^(S^ V^(^A \ <Pb ^ ) , where the cor- 
relation operator is defined as E{lp^^ ,(p^j^^) = A{lp^^^) 
B{ip^j^^). Since each local observable is bounded as 

|(^(</?A ))l < 1' < 1. we obtain a Bell inequal- 

ity from the expectation value of B: 



\Bn\ = |Tr[a 



{N+N' 



B]\ 



(6) 



i?(^i'\^g^)|<2. 



where the correlation function is 

E{v^a\'Pb^)= Y1 e{nc,mc)x 

{7ia+mc+nd+mo=N+N'} 

e{nd, mo) P^'^^'-'^^'^'^cmc; ndmo). (7) 

Here P^'f^A >v'~b ') {n^rac] ridmo) is the joint probability for 
the case that the outcome of Alice and that of Bob 
is the trace of the projection operator onto the states 
\nc,mc)cC and \nd,rn£))dD for the measurement setting 
Lp'^^ for Alice and ip^g for Bob. Therefore, if we can 
demonstrate the violation of the given Bell inequality ^ 
for a quantum state, (t(^+^ ) ^ then we can conclude that 
the state is entangled and, if the conditions for locality 
are met, also non-local. 
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IV. APPLYING THE BELL INEQUALITY TO 
PAIRS OF MODE ENTANGLED STATES 



Here we apply our Bell inequality (|6| to the pairs of 



states, \'ipN 



,«i2 



We will now evaluate the Bell inequality (|6]) with dif- 
ferent mode entangled states. In the following section, 
we apply the Bell inequality to pairs of states of a zero 
temperature, non-interacting Bose Einstein condensate 
of fixed number of particles. In section (IVB), we apply 



the Bell inequality to states that are useful for precision 
measurement, such as the so called NOON states and the 
'spin' squeezed states. 



A. Non-interacting Bose-Einstein condensate 

The zero temperature state of a non-interacting Bose- 
Einstein condensate of fixed particle number that is sym- 
metrically distributed between two modes is 



1 ^ 



\n,N -n). 



(8) 



Let us first consider the case of just a single 
particle {N = 1), so that the composite state is 

|V^i)»2 = / _i^(|10) -I- |01))j .The correlation function 

is Ei,^^{^'^^,^%^) = sin2((0^) - 0(^'^)/2), with the cor- 
responding Bell term, i37v=i, constructed via Eq. (|6|. 
When -Bjv=i > 2, the state distributed between Alice 
and Bob is non-local. The left hand plot in figure ^ 
shows on a violation of Bj^^i for a range of measure- 
ment settings. We will discuss the results in more detail 
in the following section. 



N 



N = 2 



N = 3 




with the state 



0, ^^A 



FIG. 2: Three plots of the Bell terms, Bm=i, Bn^2 and Bn^3, 
has a maximum violation of Bn^i = 2.41 at measurement settings: 
quantity Bn=2 has a maximum violation of Bn^2 = 2.36 at measurement settings: 



from left to right. Bell quantity Bm=i 
1.57, </)^^' = 3.93 and t/)^' = 2.36. Bell 
0^^' = 0, = 1.07, (j>^g^ = 3.68 and 



— 2.60. The final plot (on the right) is the Bell quantity Bn^3- Here the maximum violation is _B = 2.24 at the angles of 



= 0, 



= 1.00, 



4^' 



= 3.64 and 



4^' 



= 2.1 



We next consider N = 2 particles in each system, so 
that the composite state is |V'2)'*^- Here the correlation 

function is En=2{^^a\^b'') = sin'^ii<l>A " 05^2), from 
which the Bell term, Bpf^2^ can be constructed. The cen- 
tre plot in figure ^ shows the Bell term, Bn=2- While 
there is a violation of the Bell inequality, compared to 
the N ~ 1 plot the maximum violation is smaller and 
the range of measurement settings that give a violation 
has also reduced. For N > 2 particles, the correlation 
functions become more complicated. For simplicity, we 
show only the plot of the Bell term, Bjy^^, for the com- 
posite state, Itps)'^'^, which is on the right hand side of 
figure (pi). Again the magnitude of the violation and 



range of measurement settings that give a violation are 
smaller than in the cases of l^i)^"^ and |V'2)®^- Note that 
for all the Bell inequalities in this section the value of the 
Bell term depends on the relative measurement settings 
between the parties, — ipg and not on the absolute 
phases of the two states. 

We have also checked our Bell inequality with pairs 
of states, IV'at) ® \'^N'), where N ^ N' for balanced 
beamsplitters (a = /? — ^) on both sides and 

found the correlation functions are E{(p'"2' , tpg^) = for 
(N, N') = (1, 2), (1, 3), . . . , (1, 6),^ (2, 3), . . . , (2, 6), (3, 4) 
irrespective of measurement settings. For unbalanced 
beamsplitters (a 7^ f3) on both sides, in general 
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E{'£j^ ,ipg) ^ 0, but we have still found no violation 
of m over all parameters for (iV, N') = (1, 2). Note that 

E(^'a,Vb^) for {N,N') is equal to that for (iV',iV). 

B. Highly entangled states 

We now apply the Bell inequality to states which con- 
tain more entanglement compared to the state of a non- 
interacting Bose Einstein condensate. 

1. NOON states 




FIG. 3: The left hand plot shows the identical Bell terms, 
BN=2,m^o and -Bjv=4,m=i, with the state ([9|. The maximum 
violation here is -B]v=2,m=o = BM=3,mi = 2.41 at one of the 
four measurement settings: <^^' = —0.13, (f>^^^ — 0.65, cj>g^ = 
0.26 and <^^' = -0.52. The right hand plot shows the Beh 
term, _Bjv=3,m=o, whose maximum value is _Bjv=3,m=o ~ 1.71, 
hence there is no violation for all measurement settings. 

The so-called NOON states and their generalisations 
have been used, among other things, to gain drastic im- 
provements in precision measurements |30H32| over the 
standard quantum limit. Here we check our Bell inequal- 
ity Q with pairs of the following states 

\N,m) = ^{\N -m,m) + \m,N -m)). (9) 



Note that these states were also called MssM states in 
|32j standing for '^many-some + some-many" . 

We first consider the composite state, \N ^ 2, m ^ 
0)^2 ^ (1/72(120) -I- |02)))®^ and compute the corre- 
lation function to be E^=2,m=o{v'A ^ ^b^) — ^^^{'J'a ~ 
ip^^)- The individual state, |iV = 2,to = 0), was cre- 
ated in j33j via second order tunneling. The correspond- 
ing Bell term, B]\;^2,m=o, is plotted in figure ([3|. As 
with the case of Itpi)^"^, the maximum violation nere is 
BN=2.m=a — 2.41; however, the number of regions of 
violation has increased from one to four, which corre- 
sponds to the enhanced phase sensitivity that such a 
state would bring in a precision measurement. More- 
over, the Bell term BN=2,m=o is identical to i3Ar=4_m=i 
for the state \N ^ 4, m ='l)«52 = (l/^(|31) + \13)))'^^. 
Conversely, we have checked our Bell inequality for the 
states, |7V = 3, m = O)®^ ^ (l/V2(|30) + |03)))^2^ 
= 3, m = 1)«52 ^ (l/V2(|21) + |12)))®2 ^ 
4, TO = 0)®2 = (1/V^(|40) + |04)))«'2 and we have found 
that there is no violation for all measurement settings. 




FIG. 4: The squeezing parameter, Es, (left hand side) and the 
maximum violation of the Bell inequality ([6| (right hand side) 
for c < 1/V2. The state is spin squeezed when Es < 1. Our 
Bell inequality detects weakly entangled states (0 < c < 1/2) 
where the squeezing inequality pO] ) does not. 



2. Squeezed states 



Spin squeezing |341 135) is a mechanism that generates 
states that surpass the standard quantum limit in preci- 
sion measurements. Here the N particles, that are dis- 
tributed between the two spatial modes, are described by 
a fictitious J = N/2 spin 155] . The two modes would rep- 
resent the two states required to perform interferometry 
if that were our objective. 

Spin squeezing is achieved when the fluctuations in 
one angular momentum direction are reduced, while the 
coherence is preserved in at least one of the other two 
directions. We take here Sz = (l/2)(a^a — b'^b), Sy — 
(V2)(atS- 6to) and = (l/2)(a^6 -h 6^a), where a and 
b are the annihilation operators for modes, a and b. The 



squeezing between the two modes is given by 



(10) 



with Es < 1 corresponding to a spin squeezed state. 

While in experiments [37] spin squeezing is generally 
achieved via a non-linear interaction with on order of 
10^ particles, here we consider a toy example, which will 
allow comparisons to the other results in this section. 
We apply our Bell inequality to the composite state, 
\Mc))^^, where 



\Mc)) = c|20) + Vl- 2c2|ll) + c|02). (11) 
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FIG. 5: The Bell terms (|6| for the state ([TI| for c = 0.6, c = 0.65 and c = 0.7 as a function of two of the measurement settings 
(the other two are fixed to maximise the Bell term). The maximum violation increases with increasing squeezing, namely 
Bn^2{c = 0.6) = 2.394, Bm=2{c = 0.65) = 2.405 and Bn=2{c = 0.7) = 2.413. The number of regions of violation grows with 
increasing squeezing, which is an indicator that the standard quantum limit in a precision measurement would be overcome. 



The parameter, c, controls the amount of squeezing be- 
tween the two modes; when c — 1/2 the state is that 
of a non-interacting Bose- Einstein condensate, |'02), and 
when c = l/^/2 the NOON state, |7V = 2,m = 0), from 
earlier in this section is reached. The squeezing parame- 
ter. Eg, is plotted on the left hand side of figure Q for 
different values of c. 

Figure ^ shows the Bell term ^ for three different 
values of c of increasing squeezing, namely c — 0.6, c — 
0.65 and c — 0.7. For c = 0.6, the Bell term still be- 
haves in a similar manner to the N = 2 non-interacting 
case from figure ([2|; there is one region of violation, but 
the squeezing has increased the maximum violation from 
Bjv=2 = 2.36 to Bm=2{c = 0.6) = 2.394. The land- 
scape of the Bell term, however, changes considerably 
as the squeezing gets larger still. For c = 0.65, there 
are two clear regions of violation indicating the ampli- 
tude of the 1 11) subspace is decreasing and the phase- 
enhancing {|20), 1 02)} subspace is playing a more sig- 
nificant role. The maximum violation for c = 0.65 is 
Bjv=2(c = 0.65) = 2.405. The final plot on the right- 
hand side of figure Q shows the Bell term for a state 
which predominantly consists of |20) -|- |02), with a very 
small amount of 111) included. Here we see four clear 



domains of violation and a maximum violation that is 
identical to the state |iV = 2, to = 0). 

C. Weakly entangled states 



We can also apply our Bell inequality to weakly en- 



tangled states by taking the state (111 and allowing the 
parameter, c, to be set below the value for the non- 
interacting Bose-Einstein condensate case, i.e. c will be 
less than 1/2. By applying the standard von Neumann 



entropy to ( 11 1 one can detemine that the state, |V'2(c)), 
becomes separable only when c = 0. Here the state is 
1 11), which would correspond to the Mott regime of the 
Bose-Hubbard model. 

Figure ^ shows the Bell Term ^ for the state 
|V'2(c))'*^ for values of c where the entanglement is less 
than in the non-interacting BEC case. For all values of c 
there is just one region of violation. As the entanglement 
in the state decreases, so does the maximum violation as 
can be seen more clearly in figure Q . The Bell terms also 
become increasingly flattened with decreasing entangle- 
ment. When c = the Bell term is flat and there is no 
violation of the Bell inequality. 



ci = 0.1 



c-i = 0.2 




FIG. 6: (Color online) The Bell terms Q for the state (111 for c — 0.1, c = 0.2, c = 0.3 and c — 0.4 as a function of two 
of the measurement settings. The entanglement of the state increases from left to right as does the range and the size of the 
violation of the Bell inequality. Specifically the maximum violations for the four plots are B{c = 0.1) = 2.032, B{c = 0.2) = 
2.116, B{c = 0.3) = 2.220 and B{c = 0.4) = 2.307. 
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V. DISCUSSION 

In the previous section, we applied the Bell inequal- 
ity (|6| to various mode entangled states of fixed parti- 
cle number. For the case of the non-interacting Bose 
Einstein condensate in section (IV A), there was a de- 



crease in both the size of the violation and in the range of 
measurement parameters over which a violation occurred 
with an increasing number of particles. Restrictions on 
the measurement space due to the particle number su- 
perselection rule prevent the maximum violation (2-\/2) 
of ^ occurring for any particle number, N, over any 

set of parameters, {(p^^\ ip^^\ (pg\ ipg^ } . We can explain 
this in the following way. 

Let us take the state, IV-i)®^ = (75 (|10) + |01)))®^ 
If on each run of our Bell test one particle were always 
guaranteed to end up with each party, then the max- 
imum violation would be obtainable with the correct 
choice of measurement parameters. One can see this by 
applying the Bell inequality ^ to the renormalised state 
found by projecting iV'i)'*^ onto the subspace spanned 

by {|1010)aAbB, |1001)aAbB, |0110)aAhB, |0101)aA6i3}. 

However, in practice there is always a finite probability 
for one party to detect both particles on their side of 
the system, which implements a local particle number 
measurement on each of the four modes a. A, b and B 
and will not distinguish between the quantum and clas- 
sical correlations (see section ([ll]) for a brief discussion 
of this point). It is this mixing of outcomes from the two 
different measurement spaces that stops the Bell terms, 
-Bjv, from reaching their maximum value. 

On the other hand, we know from Gisin and Peres 
[39j that for a spin singlet state of any size, s, one can 
always find measurement settings that give the maximum 
violation of 2-^/2 to an inequality identical to The 
measurement operators used in [39] , that gave rise to the 
maximum violation, are block diagonal; with each block 
consisting of a 2 x 2 rotation matrix, Ry (a) = cos a + 
ax sin a, where and ax are the usual Pauli operators. 
This contrasts with our measurement operators ^ that 
are block diagonal in the n x n subspaces of constant n 
particles. Since, due to the superselection rule, we cannot 
rotate our measurement operators with transformations 
that mix the subspaces of different particle number, it is 
impossible to reach the measurement space used in [39] 
and hence it is also impossible to obtain the maximum 
violation of the Bell inequality. 

Indeed, for a class of quantum optical down-conversion 
Bell tests, Popescu et al. [40j analyzed the CHSH in- 
equality when the measurement space also included out- 
comes that were unfavourableand they found that maxi- 
mum violation of the inequality was 1 + 72 w 2.41. This 
is identical to the maximum violation we obtained in fig- 
ure (|2]) for TV = 1 and for the = 2, m = 0) state. 

The decrease in the the range of violating measure- 
ment parameters with an increasing number of particles 
is also due to the restrictions on the measurement space. 



This can be compared to an early result by Mermin [35] , 
who considered a Bell inequality for pairs of spin s parti- 
cles. There Mermin created a restriction on his measure- 
ment space by considering only Stern-Gerlach type de- 
vices whose operation depends solely on the orientation 
of the quantization axis and thus cannot make projec- 
tions onto arbitrary states of the subsystems. In agree- 
ment with the results in section ( IV A I , Mermin found 



that as the size of the spin increased, the range of angles 
for which the contradiction arose decreased. Speaking 
somewhat loosely, this restriction is similar to the fact 
that here when one party receives M particles, the re- 
maining 2A^ — M particles are always detected by the 
other party irrespective of the measurement setting cho- 
sen by each party. Certain combinations of measurement 
outcomes are just impossible. 



In fact, Wiseman and Vaccaro suggested in [T^ that to 
correctly determine the entanglement for superselection 
rule restricted states one should first project such states 
into the subspace of fixed particle number, calculate the 
von Neumann entropy for resulting renormalised states 
and then take their average. If one does this for pairs 
of states, IV'Af)'*^, the amount of entanglement peaks for 
N — 2 and then goes to zero by iV = 9. Since we are 
also measuring in the subspaces of fixed particle number, 
we should likewise not expect to detect any entangle- 
ment for higher numbers of particles. In order to see a 
maximum violation, one would need to make arbitrary 
measurements on the modes by coupling, for instance, to 
a Bose-Einstein condensate reservoir as in [5TJ [5S] . 



Our Bell inequality only shows violations for the pairs 
of NOON states with = 2, m = and A^ = 4, to = 1. 
This is due to the measurements in our inequality being 
linear in particle number and including no higher order 
correlation functions. As the basic group [41' that forms 
the correlations in the NOON states increases, i.e. for A^ — 
2to > 2, one would need second order observables (and 
higher) to detect the correlations. In contrast to the non- 
interacting case, the entanglement of pairs of the NOON 
states as measured by Wiseman and Vaccaro's projected 
von Neumann entropy [T^ remains constant for all A^, so 
that in principle the nonlocality of these states should be 
detectable within a different scheme. 



We also checked the CGLMP inequality [H] with the 

joint probabilities, P^'^a ''^'B'){n^rnc\ ndniu), for the non- 
interacting Bose-Einstein condensate and found no vio- 
lation. This suggests that our Bell inequality is particu- 
larly suited for detecting entanglement of pairs of states 
restricted by superselection rules. Moreover, our Bell in- 
equality is able to detect the mode entanglement of even 
weakly entangled states of two modes, A and B i.e. in 
state \tpN), when the well known spin squeezing inequal- 
ities [351 137] do not. 
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VI. IMPLEMENTATION WITH MASSIVE 
PARTICLES 

Finally, we discuss how to test this Bell inequality in 
realistic conditions. For massive particles one can cre- 
ate the mode entangled state, IV'w), by cooling TV bosons 
into the ground state of a double well potential [Ml [45] . 
Wc require an identical pair of such systems for our Bell 
test. The double wells would be positioned so that to- 
gether they form a square like shape (see figure ([T]) for 
a rough indication of the set-up). The potential barriers 
between each of the wells would initially be high, while 
maintaining the coherence of the particles. 

To implement the beamsplitting operation each party 
lowers the potential barrier between for their wells for a 
desired time depending on the beamsplitter coefficients, 
a and /3, generating an exchange of particles between 
modes, a and A |46j . Similar beamsplitter networks were 
used in [57] to detect multipartite entanglement between 
bosonic particles (as opposed to between bosonic modes, 
as in this paper), and subsequent work [48' showed that 
high precision beamsplitters for ultracold bosonic atoms 
can be realized based on current experimental technol- 
ogy. In our test, the different measurement settings are 
controlled by changing the relative phase between the 
two modes, a and A, for Alice and b and B for Bob, by 
each party locally changing the bias of one mode (poten- 
tial well) relative to the other, by applying, for instance, a 
dispersive laser pulse for a desired time |46) . The number 
of particles in each of the wells is then counted. While it 
is at the moment difficult to resolve different numbers of 
massive particles, steps in this direction have been made 
P5] . Each party would randomly choose a different mea- 
surement setting on each run of the test and the resulting 
value of the corresponding Bell term would be generated 
statistically over many runs. Note that for massive par- 
ticles it is unlikely that the measurements here would be 
performed at a speed faster than any communication be- 
tween the modes, so that locality would not be assured, 
as is actually the case for all recent Bell tests with massive 
particles (note that a proposal was recently put forward 
for a loophole free Bell test using massive particles via 
entanglement swapping f43j ) . 

We note that one could also test this Bell inequality 
with mode entangled states of photons, which would al- 
low the locality loophole to be closed. To generate the 
entangled state, IV'at), one would send two lots of N pho- 



tons through two 50:50 beamsplitters, one output of each 
would be sent to Alice and the other output of each to 
Bob's side of the experiment. These outputs would in 
turn then be passed through another beamsplitter, one 
for Alice and one for Bob, each with the desired reflec- 
tivity. To switch between the measurement settings each 
party would pass one of their modes through a different 
phase-plate prior to the final beamsplitter. Photon num- 
ber would then be measured by each party in the output 
ports of the second set of beamsplitters and the magni- 
tude of the Bell term generated over many runs of the 
test. 



VII. CONCLUSIONS 

In this paper, we have addressed the problem of per- 
forming a Bell inequality on quantum states that are re- 
stricted by the particle number superselection rule. In 
particular, we have focused on states with a fixed number 
of massive particles, so that the particle number super- 
selection rule is in effect. We derived a Bell inequality 
that allows to by-pass the superseletion rule - in order to 
perform measurements other than local particle number 
measurements, two copies of the states are used. We test 
the Bell inequality with different mode entangled states 
and find that for a non-intereacting Bose-Einstein con- 
densate, while the violation is not maximal, we detect 
some entangled states that the CGLMP Bell inequality 
cannot. Moreover, the Bell inequality presented here can 
detect not only spin squeezed mode entanglement, but 
also relatively weak mode entanglement. Our Bell in- 
equality can be implemented with current technology. 
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